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INTRODUCTION 

a 

S In the presence of strong interactions, spin and charge of electrons behave dif- 

ferently. In the half-filled case of the Hubbard model with large on-site repulsion, for 
example, charge excitation has a gap but spin excitations remain gapless. The dynamics 
in the large energy scale can be understood in terms of perturbation theory. However 
study of the infrared behavior requires more sophisticated approach like the renormal- 
ization group or exact solution. Concept of the fixed point turns out to be vital in 
understanding the dynamics in the infrared limit. Although starting models have an ex- 
treme variety, the number of different fixed points are very few. Therefore classification 
of models according to their fixed points are useful. 

The dynamics characteristic of a fixed point is in general restricted only to a narrow 
energy range, and different behavior sets in as the excitation energy is increased. Even 
very near the fixed point, presence of marginally irrelevant operators causes complication 
in dynamics and in asymptotic behavior of correlation functions. This is also the case 
for the Hubbard model and the short-range t-J model in one dimension. Thus it is 
illuminating to investigate models with simpler but nontrivial fixed-point behaviors. 
Recently a new class of interacting fermion models have been found where the fixed- 
point behavior persists to all excitation energies. 1-6 A fundamental role in this surprising 
phenomena is played by the supersymmetry which is associated with transformation of 
a fermion into a boson and vice versa. 

In this paper we discuss spin and charge dynamics in a SU(z/) generalization of 
the one-dimensional t-J model with long-range interactions. In this model dynamics 
of spin and charge are independent of each other for all energies. Moreover there is a 
discontinuity in the momentum distribution in the ground state, in contrast to the power- 
law singularity which is ubiquitous in one dimension. We have ascribed these properties 
to realization of the fixed point of free Luttinger liquids for spin and charge. 1-3 If one 



takes the high-density limit and freezes the charge degrees of freedom, the model is 
reduced to the long-range exchange model proposed by Haldane and Shastry. 7,8 

Complete integrability of a group of long-range models has been proven by finding 
infinite number of conserved quantities. 9 " 13 Operator algebra has been very powerful 
in the proof, but has been less so in deriving wave functions. As a complementary, we 
take in this paper a more primitive approach which can derive wave functions explicitly. 
In our representation one can see clearly how the independent charge- and spin-current 
excitations influence the energy of the system. Our formalism makes maximum use of 
rotational invariance of the model as well as that of the ground-state wave function. 
The technique is rather general and can be applied to a class of related models. 14,15 

We briefly discuss also a supersymmetric continuum model 16 which constitutes a 
fixed point of free non-relativistic SU(z/) fermions. The unique feature of the latter 
model is that it is solvable in arbitrary dimensions. The model has a close relation to 
the sigma model. 17 



SU(i/,l) SUPERSYMMETRIC t-J MODEL 

The t-J model is a standard model for studying strong correlation of electrons in 
copper oxides and other related systems. In this paper we generalize up and down 
components of spins to v kinds of internal degrees of freedom, which we still call spin. 
In this SU(z/) version, the t-J model is defined by 

n = Ha E flfi* + Impv - rnn^yp, (i) 

where fj a represents the annihilation operator of an electron at site j and spin er(= 
1, . . . , v), V is a projection operator to exclude double occupation of any site, and Pij 
denotes permutation operator of spins. 

As the one-dimensional lattice we consider a ring which consists of Nl (even) 
lattice points Xi with unit spacing and impose a condition tij = Jij. In order to show 
the presence of supersymmetry in the case of = J^, it is convenient to introduce a 
hard-core electron at site i as a composite particle of a fermion with a creation operator 
cj and a boson with b ai . We represent the creation operator of a hard-core electron as 
X[° = 6^c|. The spin flipping operator from (3 to a is defined by X^ 13 = $ ai bpi. We 
further define X®° = Cjc|, which is the projection operator onto the fermionic vacuum. 
A general X-operator is represented by X ab where a and b denote either or one of a. 

The commutation rules of these X-operators at a site obeys those of a Lie superal- 
gebra associated with the supergroup SU(z/,l). 18 Hence the X-operators generate, like 
the spin operators in the case of SU(2), "superrotation" which mixes cj and b ai . Let us 
introduce a (u+ l)-dimensional vector operator ^ whose conjugate has components 

*! = (&L,---,^,/4), (2) 

where h\ = Ci creates the vacant state. The superrotation leaves the constraint I = 
^j^i = 1 invariant. Thus the supersymmetry at the single-site level is merely a rep- 
resentation of the hard-core constraint of fermions. However supersymmetry imposes a 
nontrivial symmetry for intersite interactions. We note that an intersite operator ^j^j 
is also invariant under the global [y + l)-dimensional superrotation which is common 
to all sites. Thus if the intersite interaction in a model is expressible in terms of 
there remains a global supersymmetry. 



Let us introduce a graded permutation operator Py. It is defined by 

4- = : (tfj^) : = 6) (3) 

where : • • • : indicates the normal ordering of component operators. If the ordering 
makes exchange of two fermion operators, the sign is reversed. This operation leads to 
the second equality where 0& = 1 if 6 = 0, and 8b = — 1 otherwise. The presence of 
the factor 8b is characteristic of the graded permutation. We note that is invariant 
under the superrotation as apparent from the first equality. 
The t-J model is rewritten as 

H = J2J ij (P ij -l+2Xf ). (4) 

From eq.(P, it is clear that the model without the Xf° term has supersymmetry. 19 This 
chemical potential term breaks the SU(z/,l) supersymmetry just as the magnetic field 
breaks the SU(2) symmetry of the Heisenberg model. Once the model without the Xf° 
term is solved, inclusion of the term is easily done. We note that the discussion of this 
section is general as long as £y = Jij. 



PERMUTATION PROPERTIES OF SU(i/) ELECTRONS 

The ground-state wave functions for a group of long-range SU(z/) fermion models 
have been shown to be of Jastrow form. 6 ' 15 The derivation uses elementary but compli- 
cated algebra. The complexity arises when one deals with permutation of spins. We 
note that the permutation property is in fact the same as that of free SU(V) electrons. 
By using this fact, one can bypass most of the complicated algebra in deriving the 
ground-state wave function and the energy of the SU(V,1) t-J model. 

Let us consider free N e electrons in the continuum ring with unit radius. At the 
ground state the electrons fill the momentum states from zero to the Fermi wave number 
for each spin. This leads to a spin singlet wave function ^f{8, cr} which is given by the 
product of Slater determinants. Here {8, cr} denote the set of spatial and spin coordinates 
with < 8i < 2tt for each particle. In one dimension, the Slater determinant is reduced 
to the Vandermonde determinant and ^f{8, cr} is given by 

if! F {6, a} = exp(iS'), 

1 7T 

S = XK^i, <Tj) ln ( sil1 o%) + TTSgnfo - <Tj)]. 
i<j Z Z 

This wave function has a nontrivial permutation property. To derive this we apply 
the momentum operator pi = —id/ddi to ^f{8, cr}- Then we obtain 

PiV F {0, °} = -iJ2 5 (^> Oi)cy*f-{0, cr}, (5) 

where Cjj = 2 _1 cot(8ij/2). The singlet nature of ^p{8,cr} makes it invariant against 
any SU(i>) rotation R, and obviously RpiR' 1 = Pi. Rotational average of eq.fl5|) leads to 
replacement of 5(<Ji, <jj) by /Ty where Ty is a projection of operator to spin symmetric 
states and the prefactor / = 2/(u+ 1) corresponds to the fraction of equal-spin pairs to 



the total number of symmetric spin states. In terms of the spin permutation operator 
Pij, Tij is given by = (P^ + l)/2. We then introduce an operator 

Oi Pi i ^ ] fCijTijj 

which annihilates the ground state for each i. Consequently we obtain Y^i a l a i^F = 0. 
After some algebra this relation leads to the following identity: 



/£ 'CijCu-Tij^p 



i(l-/)E^sin- 2 (i%) + ^ 



*F, (6) 



E v = -—N c {vNl + 3N a + u- 3), 

where N a is the number of electrons for each spin. The sum with prime runs over 
coordinate indices different from each other. 

This identity relates the three-body term in the left-hand side to the two-body term 
plus the constant. Two remarkable features of this identity are worth noting: First, the 
identity is valid also for the Gutzwiller-projected wave function in the one- dimensional 
lattice, since restriction of d{ on the lattice together with hard-core constraint has no 
effect on the permutation property. Secondly, the identity is also valid for any wave 
function which is a product of ^fp with a spin-independent function. 

As an analogue of the ring system for SU(z/) free fermions, one may consider a 
linear system where a harmonic potential prevents fermions from going to infinity, the 
wave function of the SU(z/) free fermion ground state $f{x, a} is easily constructed. In 
this case we obtain the analogous identity which reads 



(Xi - Xj)(Xi - Xt) 



Note that the identity is independent of strength of the harmonic potential, and there 
is no constant term. This identity is useful for exact solution of the SU(V) Calogero 
model, 14 and simplifies the treatment by Vacek et al. 15 



WAVE FUNCTION AND ENERGY OF THE GROUND STATE 

From this point on we restrict ourselves to the case where the interaction has the 
long-range form: = Jij = tD(xi — Xj)~ 2 with D{x) = (N l /ti) sm[iTx/N L } and assume 
t > 0. For this SU(z/,l) t-J model Ha and Haldane have given a nice and detailed 
account of their calculation 6 which is rather involved. We give in the following a different 
treatment which we believe is simpler. The formal proof that the wave function indeed 
is the ground state is still missing, but there are strong evidences in favor of that. 1-6 

The major role in our treatment is played by the identity eq.(|5|) for the Gutzwiller 
wave function. In order to treat the transfer term together with the hard-core constraint 
it is convenient to take a fully polarized state, say in the direction of u, as the reference 
state \Nlv). Then there emerge fermionic holes and hard-core bosons with spin a (= 
1, . . . , v — 1) as particles. The SU(z/) Gutzwiller state \G) is represented by 

\g)= e w.Min n /£/* n /^>- 

{v} c=i iefav} je{y} 



Here {x a } denotes the set of coordinates for N a electrons with spin a and {y} does 
that of Q holes. Thus we have Nl = vN a + Q in the singlet state. In order to remove 
the degeneracy we choose N a odd. The amplitude ^a({x a }, {y}) is given, apart from a 
normalization factor, by 

Vg({x<t},{v}) = ew[-in(J2 x i° + Y.ye)]nn D ( x iv - x j°) 2 

ia f. a i>j 

x \{\{D(x, ia -x jfi )\{Y{D{x ia -y,)\{D(y,-y m ).. (7) 

i,j a i,£ l>m 

In this representation it easy to apply the transfer operator of electrons with spin 
v, since it is translated into the motion of holes without any effect on Xi a . On the other 
hand, transfer of electrons with a different spin has the same effect as that of spin v 
in the singlet state. 1 Thus the whole effect of the transfer is just v times of the hole 
transfer in eq.(^). Equivalently, we can take rotational average of the wave function after 
calculating the hole transfer, and multiply the result by v. As a result of the transfer, 
we have three-body terms multiplying ~^q. At this state we take the vacant state as 
a new reference state to make use of the SU(z/) invariance. Then the three-body term 
turns out to be of the same form as the left-hand side of eq.(|6|). By using the identity 
ec L-© we obtain the two-body term which is just minus of the interaction term in eq.(|I]) 
and the constant. The ground-state energy E g corresponds to minus of this constant 
and is calculated as 

4 = ~{nl - 3n 2 + 2vn e ) + -L[(„ + 2)n e - 3i/], 

7T t 01/ bl\L 

with n e = N e /Nl = uN^/Nl- This results agrees with that of ref.6 and reduces to that 
of ref.l for v = 2. The charge susceptibility Xc is calculated from the second derivative 
of Eg as Xc = 7r 2 t(l — n e )/u. 



SPIN AND CHARGE CURRENTS IN THE FIRST QUANTIZATION 

It has been shown that long-range the supersymmetric t-J model has completely 
separated spin and charge excitations. 1-6 In this section we clarify how this separation 
is related to the particular form of the model. The single-particle spectrum would be 
given, if there were no two-body interactions, by Fourier transform of or tij. This is 
calculated as 

J(q) = £ Jij exp(-zg^) = (q- ^) 2 - l(iV 2 + 2), (8) 

i 

where 9ij = Q{ — 9j = 2n(xi — Xj)/NL and q is an integer with < q < Nl- The unit of 
energy is so chosen that J — t — (N L /n) 2 /2. 

A remarkable feature of this spectrum is that J(q) is a quadratic function of q. If 
one shifts the origin of q to the edge Nl/2 of the Brillouin zone, the spectrum is nearly 
the same as that of free particles in continuum space, except for the presence 20 of the 
cut off in q. Using this similarity we can also work conveniently with the first-quantized 
representation with the completely polarized state in the direction of a = v as the 
reference state. The complete set for the system consists of the product of one-body 
states inside the Brillouin zone with proper symmetrization. In terms of Zi = exp(i^), 
where 9i denotes a coordinate with either a spin other than v or a hole, the one-body 
states are spanned by monomials z\ with < k < Nl. We now shift the origin of the 



wave number to the edge N^/2 of the Brillouin zone. Within the many-body Hilbert 
space defined in this way, we may replace q — Nl/2 in eq.(|S]) by —id/ 89. Then the 
Hamiltonian is given by 

M 1 1 

« = + 4 E sm- 2 (-%)(l + P<j) + £m, (9) 
i=i i^j 

where = —id/dOi and M = iVx — N a is the number of holes plus bosons with spins 
different from v. The graded permutation operator acts now in the space of SU(z/ — 1) 
spins and holes. The constant Em appears as a result of transformation. 

By taking the hole picture in the first quantization, the minus sign in the hopping 
combines nicely with the plus sign in in leading to eq.(^) where pj is common to 
holes and bosons. Equation can be regarded as the SU(z/ — 1, 1) generalization of the 
Sutherland model in continuum space. Thus all eigenfunctions of the SU(V,1) t- J model 
can be mapped into those of the SU(z/ — 1,1) Sutherland model. However the reverse is 
not true because of restriction of wave functions to the Brillouin zone in the t-J model. 

We introduce a shift k a in the momentum distribution for each spin in the original 
lattice system. In the presence of small currents for each spin component a or j3, the 
energy of the system increases by the amount 

SE= — — 22 v a pJ a Jp, 



2N L 



al3 



where a current component J a = 2k a is an even integer. The velocity matrix v a p is con- 
strained to the two-parameter form v a p = S a pu + v. Here the off-diagonal component is 
independent of spin indices because of the SU(z/) invariance. By making linear transfor- 
mation of J a , one can diagonalize the velocity matrix. The first eigenvalue corresponds 
to the charge velocity given by u + v , and the remaining v — 1 ones are spin velocities 
which are all degenerate and are given by v. Thus if one can calculate SE for two sets 
of {J^}, one derives spin and charge velocities. 

In the first quantization representation, the operator to generate the current be- 
comes diagonal in coordinates and is given by 

<P{k}{9} = exp li KYl ^ ~ i J2 k ^i 

where di a refers to a coordinate of a particle with spin a, while 9j refers to all particles. 
The minus sign with k v comes from the nature of fully polarized reference state. 

There are two particular kinds of currents for which SE is easy to calculate. The 
first kind is k a = for o ^ v. Then the generator boosts all particles in the first 
quantization by the equal momentum — k u , and commutes with Pij. Then we obtain 

M M 

H<P {k} V G = *gE Pi<f>{k} + 2 E W«)ft*G + <P{k } n^ G , (10) 

1=1 1=1 

where the second term on the right hand side is zero. The increment of energy is thus 
given by 

SE = {N L - N a )kl. 

The second kind of currents corresponds to the charge current, and is given by k a = k 
for all a. In this case the exponent in <f){k}{9} cancels each other except for that of holes. 
At this point we notice that for the singlet state the transfer of electrons with spin v is 



in fact equivalent to transfer of any kind of spin. Then for the singlet state the effective 
Hamiltonian is given by 

i=l cr&v) i i^j 

where is the SU(v — 1) permutation operator. Note that the hole degrees of freedom 
appear in the first term only. We emphasize, however, that the hard-core repulsion 
between holes and spin bosons must be taken into account in the solution of eq. ( pi]) be- 
cause the effective Hamiltonian is valid only for the lattice model. With this restriction, 
we can still proceed in the same way as in eq. (|T0|) and obtain 

5E = uk 2 Q. 

We have thus two independent equations which are enough to determine the spin 
and charge velocities v s and v c . In the original energy unit they are calculated as 

v s = v = 7rt, v c = u + v = 7r£(l — n e ). 

This result is consistent with refs.l and 6. 



DISCUSSIONS 

In the present derivation of current excitation spectra, it is clear that the quadratic 
dependence of SE on currents is due the parabolic spectrum of the model. This depen- 
dence is not restricted to infinitesimal currents as long as the momentum distribution is 
within certain limit. 1-6 The parabolic spectrum is also responsible for the spin-charge 
separation for all energies. Our treatment here is restricted to the singlet ground state 
from which currents are generated. It should be interesting to introduce a generator to 
polarize the ground state, and to study the increment of energy. The generator may be 
related to the Yangian discussed in ref.22. 

In the momentum space the t-J model is written as 

N L Af 

H = - E( fc - -y) 2 Y, xab ( k ) xba (- k )^ + ^o, (12) 

k=l a,b 

where X ab (k) is the Fourier transform of Xf' and Eq is a constant dependent on the 
number of vacant sites. As we have considered polynomial wave functions of Z{ = 
exp(i8i) in the first quantization, we can also restrict operators to polynomials of X^. 
Under this restriction (q — N L /2) in eq.(B) is replaced by —id/ 06 which acts on X- 
operators. Then the model can also be written as 

h -~ Nl L 2^\ b ^r-^r db+E "- (13) 

Inspired by this form, we have recently introduced a new supersymmetric model defined 
in a d- dimensional continuum space. 16 The model is given by 

H = }igJ dr J2 VX a \r) ■ VX ba {r)6 b + C, (14) 

a,b 

where V = 8/ dr. In this continuum model one needs to perform renormalization to 
remove divergences. The constants g and C are determined so as to give zero energy 



for the vacuum and finite energy for a finite number of hard-core fermions. Note that 



In ref.16 we have shown that after renormalization all wave functions and eigen- 
values are the same as those of //-component free fermions. The presence of supersym- 
metry is essential in this surprising cancellation of the hard-core constraint against the 
exchange-type interaction. However, there is no spin-charge separation in this contin- 
uum model. 

The author would acknowledge useful discussions with N. Kawakami, H. Yokoyama 
and J. Zittartz. 
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